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O ! 

^ We calculate the persistent currents in S — -F|_F|[|] — S structures (S denotes the superconductor 

■ in a closed ring geometry, F the ferromagnet in a two dimensional geometry, and the arrows the 

alignment of exchange field) as a function of the applied magnetic flux, self consistently by solving 
O , the Bogoliubov-de Gennes equations within the two dimensional Hubbard model. The local current 

^) • shows a sign change i.e. a to tt transition in the parallel alignment of the magnetizations and it 

does not change sign in antiparallel alignment. 



o 

O ■ I. INTRODUCTION 

U 

Presently the study of tt states is a subject of intensive theoretical and experimental investigation. The tt states 
^ ' have been predicted theoretically long time ago in junctions containing magnetic impurities. They appear as a sign 
change of the Josephson critical current in SFS structures (S denotes the superconductor and F the ferromagnet). 
C3 . The crossover from the to tt state has been observed recently in the critical current versus temperature in SFS 
junctions [1], and in the critical current versus ferromagnetic layer thickness in SIFS junctions [2]. More recently 
the concept of tt state was used to explain the observed shift in the critical versus magnetic field pattern by half 
the flux quantum in the ferromagnetic — tt SQUID [3], and the decaying density of states oscillations in s-wave 
Q superconductor ferromagnet hybrid structures [4]. Similar effects have been observed in d-wave [5,6] superconductor 
O ferromagnet hybrid structures. Theoretical explanation has been given in the framework of the quasiclassical theory 
for s-wave [7] and d-wave case [8] . 
I Alternatively tt states appear in Josephson junctions involving d-wave superconductor due to the sign change of 
^ ' the superconducting order parameter [9]. For example the to tt transition occurs in the critical current versus 
temperature in Josephson junctions with d-wave symmetry and depends on the orientation of the lobes of the d- 
wave order parameter with respect to the interface. The Josephson effect has been studied in junctions between 
unconventional superconductors across different types of magnetic barriers [10]. 

In superconductor - ferromagnet multiterminal hybrid structures the proximity effect can be controlled by the 
' alignment of magnetizations in ferromagnetic electrodes connected to superconductor [11-13] or by the magnetic flux 
I in an Aharonov-Bohm loop connected to superconductor [14]. In two superconductor - ferromanget bilayers separated 
by thin insulating film, the enhancement of the Josephson critical current when the orientations of the magnetic 
moments in the ferromagnetic layers are antiparallel has been predicted [15]. Since then several S-FF-S, SF-SF 
structures have been examined in the limit of thin F layer and in the diffusive limit by solving the Usadel equations 
• [16,17]. They found that in the case of antiparallel alignment of magnetizations in the ferromagnetic electrodes the 
P5 ' critical current enhances with the exchange field, while in the parallel alignment the junction exhibits a transition to 
O • the TT state [16]. Also non-sinusoidal current phase relation in S-F-S, SF-c-SF junctions has been predicted [17]. The 

same problem has been treated in the ballistic regime [18]. 
^ In this paper our goal is to explore several new aspects related to the control of the Josephson effect in nanos- 

'kp ' tructures. We study two dimensional F-^F^^^ embedded in a superconducting wire. The method is based on exact 
. diagonalizations of the Bogoliubov-de Gennes equations associated to the mean field solution of an extended Hubbard 
^ model. The basic quantity which we calculate is the local current as a function of several relevant parameters: the 
distance from the surface, the magnetic field that penetrates the superconducting loop and exchange field. 

We find that the local current is periodic function of the fiux and is controlled by the orientation of magnetization 
in the ferromagnetic electrodes. The local current shows a sign change i.e. a to tt transition in the parallel alignment 
of the magnetizations and it does not change sign in antiparallel alignment. Our predictions from the simulations of 
this model are of interest in view of future experiments on nanoscale Josephson junctions. 

The article is organized as follows. In Sec. II we develop the model and discuss the formalism. In Sec. Ill we discuss 
metallic or ferromagnetic structures embedded inside a superconducting loop. Finally summary and discussions are 
presented in the last section. 
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II. BDG EQUATIONS WITHIN THE HUBBARD MODEL 



The Hamiltonian for the Hubbard lattice model is 

+ Vo^ni-(nii, (1) 



where i, j are sites indices and the angle brackets indicate that the hopping is only to nearest neighbors, n^cr = cJ^Qcr 
is the electron number operator in site i, is the chemical potential which is set to zero, hia- = —haz, is the exchange 
field in the ferromagnetic region and = ±1 is the eigenvalue of the z component of the Pauli matrix. Vq is the 
on site interaction strength which gives rise to superconductivity. Within the mean field approximation Eq. (1) is 
reduced to the Bogoliubov deGennes equations [19]: 

^ AW u„t(ri) \ _ / Uniin) 



such that 

iunain) = -t^^Unain + 6) + IJ,Una{n) + htazUnain), (4) 



Au„<^(ri) = Ao{ri)u„a{ri), (5) 



where the pair potential is defined by 



Aoin) = Vo < c^{ri)c^{ri) > . (6) 
Equations (2,3) are subject to the self consistency requirement 

Ao(ri) = Vo{ri)F{n) = ^ MnX^^in) tanh(/Je„^,/2)) + Uni{ri)v*^in) tanh(/3e„^,/2)] , (7) 

n 

F{ri) is the pair amplitude. We solve the above equations self consistently. The numerical procedure has been 
described elsewhere [12,13,20-22]. 

The local current between sites i,j is given by 

hj = f (* < cjcj > -t* < c]ci >) (8) 



III. RESULTS 

We study the quasiparticle properties of a two dimensional plane to which a superconducting wire is attached at 
two points and a magnetic flux is applied through the wire. The magnetic flux creates a phase difference between sites 
a and (3 (see Fig. 1), which is proportional to the applied flux. Alternatively a current through the wire can create 
this phase difference. We demonstrate in this section that the magnetic flux through the wire as seen in Fig. 1 can 
be used to control the proximity effect in this hybrid structure. The magnetic flux through the loop is modeled as a 
factor e'-^ where / = 27r$/<&o = '^''^4' in the hopping integral. $o = 7 is the normal flux quantum. In the calculation 
we used a small cluster of 8 x 8 sites to model the 2DEG and also open boundary conditions, while a chain of 20 
sites models the superconducting wire. The hybrid normal-superconducting ring has been discussed in Refs. [23,24]. 
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They found that the persistent currents depend on the length of the N segment compared to the S. For large length of 
the superconducting region the quasiparticle wave function is localized in the normal metal region and quasiparticles 
are entering the superconductor via Andreev reflection. In this case the persistent currents have period half the flux 
quantum. In the opposite limit where the length of the superconducting region is small compared to the coherence 
length the quasiparticle wave function becomes extended and the persistent currents have period equal to the flux 
quantum. In our case the quasiparticle properties are more close to the former case. We checked the case where the 
wire and the reservoir arc both normal metals. In that case we found that the persistent currents have period <I>o. 
Also the phase of the persistent currents is zero in a loop with even number of carriers and equal to tt in a loop with 
odd number of carriers. However in our case we fix the loop size and study the response of the system to the magnetic 
field and exchange field. We distinguish the following cases depending on the quality of the barrier 

A. S-N-S 

We discuss first the case where a superconducting wire is connected to a 2DEG as seen in Fig. 1. We demonstrate 

that the flux through the loop modulates in a periodic way the quasiparticle properties of the film, so that the 
proximity effect is controlled by the applied flux. The local current-magnetic flux relation is seen in Fig. 2 between 
sites aa, a(3 and jSb. It shows the characteristic sinusoidal form, with period almost $o/2. 

B. S-F-|-F-|--S 

We discuss now the case where a superconducting wire is connected to a ferromagnetic film as seen in Fig. 3. 
The objective is to investigate the effect of the exchange field in the current. The ferromagnetic atoms at sites a 
and /3 where the superconducting wire is connected to the film have the same orientation of magnetizations. Due 
to proximity effect the pair amplitude shows decaying oscillations with alternating positive and negative signs inside 
the ferromagnet away from the connection points a, /3. The local current oscillates with cj) as seen in Fig. 4. The 
Josephson effect occurs due to transfer of Cooper pairs across the interface via Andreev reflection. In the ferromagnetic 
site with spin up orientation of magnetization a and in the ballistic limit the phase difference between the electron 
and Andreev reflected hole depends on the exchange fleld i.e. Scf) = 2hx/TiVF- In the other ferromagnetic site with 
the same orientation of magnetization /? the phase shift has the same sign. So in total in the parallel alignment of 
magnetizations there is a phase shift which can be equal to tt, and the critical current is expected to change sign 
compared to the SNS case. Indeed Iai3 has opposite sign in relation to the normal metal case as seen in Fig. 4. 

C S-F-|-F|-S 

We discuss now the case whore a superconducting wire is connected to a two dimensional ferromagnetic domain 
wall as seen in Fig. 5. Wc demonstrate that the alignment of the magnetization provides additional control parameter 
for the proximity effect. The ferromagn(^tic atoms at sites a and /? where the superconducting wire is attached 
have the opposite orientation of magnetizations. The local current oscillates with (j) as seen in Fig. 6. However 
differently to the parallel alignment case, the local current lap does not change sign with the increase of the exchange 
field. These results indicate that the junction in the antiparallel alignment of magnetizations strongly resembles the 
SNS junction. We emphasize here that the ferromagnetic layers are small (one atomic layer) so there is no spatial 
variation of the pairing amplitude inside the ferromagnet. Therefore the to tt transition is explained in terms of 
phase discontinuities at the interface. For the antiparallel alignment the phase shift across the junctions is zero. 
The Josephson effect can be viewed as a transfer of Cooper pairs across the interface via Andreev reflection. In the 
ferromagnetic site with spin up orientation of magnetization a and in the ballistic limit the phase difference between 
the electron and Andreev reflected hole depends on the exchange fleld i.e. = 2hx/hvF- In the other ferromagnetic 
site with spin down orientation of magnetization (3 the phase shift is the opposite. So in total in the antiparallel 
alignment of magnetizations there is not a phase shift, and the local current is expected to be similar to the SNS case. 
For the parallel alignment the phase shifts are added and this provides the possibility for a tt state. This tt phase 
shift in the phase of the order parameter is responsible for the sign change of the current. Similar analytical results 
have been obtained using the Usadel equations [16,17]. However in that case several approximations were used like 
constant order parameter in the superconducting electrodes and strongly coupled SF interfaces. 
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IV. CONCLUSIONS 



We calculated the local current for normal and ferromagnetic planes with superconducting mirrors, in the presence 
of an external magnetic field within the extended Hubbard model, self consistently. The local current, is periodic 
function of the magnetic flux that is applied through the superconducting wire. The local current versus exchange 
field is further controlled by the alignment of the magnetizations in the ferromagnetic electrodes. Specifically it shows 
a sign change in the parallel alignment of the magnetizations but it does not change sign in the antiparallel alignment. 

Our results are influenced from the details of the lattice structure and quantitative differences are found compared to 
the results of quasiclassical theory. The results can be generalized to extended contacts. We expect only quantitative 
differences. In the present paper we used superconducting rings where the applied flux creates a phase difference in 
the interruption points. However the results can be generalized to junctions instead of rings, where the current creates 
a phase difference between the two superconductors. 
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FIG. 1. The SNS junction of the superconducting wire with the 2DEG. The loop contains 20 sites while the 2DEG is of 8 x 8 
sites. The labeling of several sites close to the interface is shown. 
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FIG. 2. The local current lij between sites i,j of the SNS junction, where ij — aa, a/3, [3h as a function of the magnetic 
flux. The exchange field is equal to ft = and the value of the superconducting gap is Vb = —3.5. 
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FIG. 3. The S-F|F|-S junction of the superconducting wire with the ferromagnetic 2DEG. At the connections sites a, 13 the 
ferromagnetic atoms have the same orientation of magnetizations. The loop contains 20 sites while the 2DEG is of 8 x 8 sites. 
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FIG. 4. The local current Uj between sites i,j of the S-F|F|-S junction, where ij = aa, af3, I3h as a function of the magnetic 
flux. The exchange held is equal to ft = 3 and the value of the superconducting gap is Vo — —3.5. 
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FIG. 5. The S-F^F^-S junction of the superconducting wire with the ferromagnet with antiparallel orientations of the 
magnetizations. The loop contains 20 sites while the 2DEG is of 8 x 8 sites. 
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FIG. 6. The local current Uj between sites i,j of the S-F|F|-S junction, where ij = aa, af), [3h as a function of the magnetic 
flux. The exchange fleld is equal to ft = 3 and the value of the superconducting gap is Vo — —3.5. 
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